Real space Dynamical Super Cell Approximation for interacting disordered systems. 
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Effective medium super-cell approximation method which is introduced for disordered systems is 
extended to a general case of interacting disordered systems. We found that the dynamical cluster 
approximation (DCA) and also the non local coherent potential approximation (NLCPA) are two 
simple case of this technique. Whole equations of this formalism derived by using the effective 
medium theory in real space. 
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Theoretical understanding of strongly correlated sys- 
tems such as high temperature superconductors, heavy 
fcrmions and magnetism require appropriate techniques 
to obtain their physical properties. Recently, by the 
many body theory techniques, it has been shown that 
both alloys and strongly correlated systems in the infi- 
nite limit dimensions are mapped to a single site which 
embedded in an effective mediumi*2iii4. This can be de- 
scribed by restriction on locality of self energy, that is, 
Yi{i,j;iu> n ) = H[uj n )§ij. This single site technique for 
strongly correlated systems and alloys systems called, 
Dynamical Mean Field Approximation (DMFA)- and 
Coherent Potential Approximation (CPA) 6 respectively. 
The single site nature of the infinite dimensional limit is 
imply that the inter sites correlations and also inter-site 
multiple scattering are negligible. But in the ordinary 
dimensions such as one, two and three, both inter sites 
correlation and inter-site multiple scattering have a sig- 
nificant contribution on the self energy, therefore not only 
the self energy is not local but also it is very sensitive with 
respect to the dimension. Recently by coarse graining of 
the self energy in the k-space, the non-local correction to 
the DMFA has been considered, this method called the 
Dynamical Cluster Approximation (DCA) 7 '®. But what 
is the correspondence of this k-space coarse graining in 
the real space was unknown. Also for a disordered system 
we have introduced a new real space technique that called 
effective medium super-cell approximation (EMACA) 9 . 
It has been shown that this method leads to the coarse 
graining of self energy and also to the average Green func- 
tion in the k-space. Also EMSCA recovers CPA, where 
number of sites in the super-cell is one (N c = 1), and is 
exact for the case of N c — * oo. Now we should ask this 
question: dose it exist a unique EMSCA, such that treats 
both interacting and randomness together?. We answer 
this question by introduce a general version of the EM- 
SCA, where in a especial case of only interacting, it is 
reduce to the DCA 7 . and for disordered systems leads to 
the old EMSCA&ia. 
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In this paper the real space effective medium super-cell 
approximation (EMSCA)2*A£ is generalized and extended 
to treat a random interacting system. We show that ne- 
glecting of the, interaction between electrons in different 
super-cells, multiple scattering by sites in different super- 



cells and also deviations of the hopping integrals, 6t° 
naturally leads to the super-cell periodicity of the self 
energy, T,(i,j;iw n ), with respect to the super-cell trans- 
lation vectors, r^. This provide us the DCA£ coarse 
graining of the self energy in k-space. We are obtain 
a closed set of equations to do calculations in a general 
EMSCA. Also whole relevant DCA equations are derived 
by the EMSCA method. 

We start our investigation on a general tight binding 
model for an interacting alloy system, which is given by, 

H = - V i- CT 'ct c / 

/ 4 y Ha^jcr 
ijcrcr' 



where c ia (ci a ) is the creation (annhilation) oprator of 
an electron with spin a on lattice site i and hi a = c^ ia Ci a 
is the number oprator. t^ are the random hopping in- 
tegrals between i and j lattice sites with spin a and a 
respectively, fi is the chemical potential and Si is the 
random on-site energy, where takes —8/2 with probabil- 
ity 1 — c for the host sites and 5/2 with probability c for 

impurity sites. U?? is a positive or negative interaction 
potential between electrons on the lattice site i and j. 

The equation of motion for electrons corresponding to 
the above Hamiltonian, Eq^ is given by, 

£((&-*i + M)M^ ~tf )^ V (ir,,V) + 

la" 

Y,U?/' G i' a ' \lT,lr,lr + ,jr) = S(r -r')%^ (2) 

la" 



where r and t are imaginary time, C" (it, jr ) is the 
random interacting single particle Green function and 

" {it, it, ir + , jr ) is the two particle Green function. 

The random hopping integrals, , can be defined in 



terms of clean system hoppings, t^ a , and the hopping 
integral deviations, 5tf a , in a such way that the hopping 
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randomness is includes just in the Stff , where 



and S(k; iuj n ) is the self energy matrix, 



,aa ,0aa , Sx(T(7 

Hi Hi t~ o^ij 



(3) 



The Dyson equation corresponding to Eq[21for the aver- 
aged Green function, G cnj (It ,jr ), iaii, 



la 

]T I drY?*' (itM)G°'° (\t\ 3 t) = 

la" J 

SiT-T)^^ (4) 

where the self energy, £ <TCT (ir, It ) , is defined by, 
J2{((^uS^+St^")G^"^(lT,jT) + 

la" 

J2u^"g{°\It,It,It + ,jt')) = 

la" 

£ [ dTTr r '\iT,lT)&"' r \lT",3r). (5) 

la" J 

The imaginary time Fourier transform of Eq0] leads to, 
]T ( + fi)5ii5 aa " - t^r" - £™" (», /; «.;„) 



/(J 



where ui n = ^(2rt + arc the Matsubura frequencies. 
Also the space Fourier transform of the EqfHJis given by, 



]T ( («„„ + M )«5 CTCT » - e--" - (k; iu n ) ) x 

a' 

G ffV (k;zu; n )=<5 CT</ , (7) 



where £ <T<T (k; iw„) is, 



S CTCT ' (k; iw„) = - (*, J'i m.)' 



(8) 



and also the clean system band structure, e£ CT is defined 
by, 



(9) 



The matrix form of Eq0in the spinor space can be writ- 
ten as, 

G(k; iuj„) = ( («*>„ + £k - £(k; «^«) ) 1 • (10) 
where the band structure matrix, ek, is 



,TT ,U 
, , / e k e k 

e k e k 



(11) 



n^n)=(f^ n \ ylit^l), (12) 

\ (k;iw„) S u (k;iw„) y 

and the average Green function matrix, G(k; iui n ), is de- 
fined by, 

G(k^)=(^^ %[^\), (13) 

also, I is a 2 x 2 unitary matrix. Hence the clean non- 
interacting Green function is given by, 

G (k;«j n )=((tw n + A«)I-6 k )~ 1 . (14) 

So, the real space correspondence of Eq[7|can be written 

as, 

G(i,j; vjj n ) = G (i,j;tu„) + 

£„, G °(h l; tu n )V(l, l;iui n )G(l',j; iw n )(15) 

where the real space clean non-interacting Green func- 
tion, G°(i,j;w n ), is, 

G°{i,j;iu n ) = I^ e -*.r«GQ(k ;i a; n ), (16) 

k 

and the average single particle Green function is defined 
by, 

G(i,j; iu) n ) = I^ e -« k - r «G(k;^ n ). (17) 

k 

EqEl can n °t be solved exactly. We extend the effec- 
tive medium super- cell approximation (EMSCA) which is 
recently introduced by us for a disordered system^, to the 
case of interacting disordered systems. The EMSCA for 
such system based on three assumptions: first, neglecting 
interaction between electrons on different super-cells, 

U?? =0, if i and j ^ same super — cell, (18) 

second, neglecting hopping integral deviations, 8tf? , 
when i and j are in the different super-cells, 

St"" = 0, if i and j £ same super — cell. (19) 

and finally, neglecting multiple impurity scattering and 
also correlations between different super- cells. These con- 
ditions imply that no correlation between super-cells, 
hence we have, 

S sc (i, j; ioj n ) — 0, if i and j same super — cell. 

This means that the self energies in each super-cell are 
independent of other super-cells and they are periodic 
with respect to the super-cell translation vectors, r^c, 

S sc (r/j + r Nc ; iu n ) = S sc (r 7J ; iu n ), (21) 
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where / and J refer to sites in a same super-cell. Also 
Eq^| leads to the super-cell periodicity for the interac- 
tion potential matrix, XJ SC (I, J), where 

TT ( T J\ ( Ujc(I,J) Ujc(I,J) \ fOO\ 

UaB ^ J) -\U^I,J) UM(I,J))> {22 > 

which is, 

U sc (r/j + r Nc ) = U sc (r/j). (23) 
The Fourier transformation of Easl21landl23limplv that 9 , 

(24) 



e -jk.rjvc _ ^ 



Thus the wave vectors of the self energy, E(k; iui n ), and 
also of the interaction potential, U(k) are restricts to the 
K„ which are given by, 



(25) 



where bi are the reciprocal lattice primitive vectors, Nci 
are the number of the lattice sites in a super-cell at aj 
direction (a; are the lattice primitive vectors) and Zj is a 
integer number. Therefore by inserting Eas l25l and 1201 in 
to Eq[S]we found that, 

£ sc (K n ;«u; n ) = -^^Te lK - r "£ sc (/, J;uj n ). (26) 



/../ 



Also by applying the EMSCA conditions, EaslT%landl2*3l 

to the following exact relation, 



(27) 



we found a similar coarse graining for the interaction po- 
tential, 

u ( K ») = ^E u -( / ' J ) elK "- r/J - ( 28 ) 



/,./ 



The inverse Fourier transformation of the self energy, 
S(K„;iw„), to the real super-cell is, 

E sc (/, J; iu n ) = -^rJ2 e~' K ~ r ' J E(K„; iw n ). (29) 



K„ 



and also for U(K„) is, 

U sc (/, J) = -l^U(K„)e- ,K - r - (30) 



K, 



where the orthogonality condition in a super-cell is given 

by, 



K„ 



(31) 



Now by inserting EqE3and[H|in to EatTTl the super-cell 
average Green function, G SC (J, J; ^w Tl ), is given by 9 , 

G sc (7, J;iiu n ) = ^£V K "' r «G(K n ;su; n ) (32) 



K„ 



where 



^E( G°(K n +h';tuj n )- 1 -i:(K n ;tuj)) \ (33) 



G(K„;iu„) 



Nc 



k' 



In order to obtain G SC (J, J; iuj n ) from Eg 1321 we should 
have S sc (/, </; iu> n ), thus we need to have another equa- 
tions to complete the self consistency loop. These equa- 
tions are obtained by applying the EMSCA to the system 
partition function, as follow. The partition function of 
the system with the Hamiltonian EqQ is given by, 



Z=(Tr e-? H } r , 



(34) 



where ( } r denotes the configurational average over ran- 
dom energies, £j. The partition function, Eg 1341 can be 



rewrite as 



Z = ( J V^V^e~ s ) r , 



(35) 



where the action S is, 



S = E j W*(W*A^-ti+%° )MT)+Sr-i, 

(36) 

and S r _i is, 



S '-i = E J dr^ ia U^ $ jgf (r)^v (r) 

ijffff 

+ l ' dT^iairfeiSijl/jja^) (37) 
ija 

+ E / ^ iCT (r)^/^ V (r), (38) 



l j a a 



in which Dty = Wid^i^d^ic and T>^ — Hidipiadtpia, 
where dtp ia = limM-»oo n^f =1 d^ CT (r m ), dif) ia - = 
limM-»oo n^f =1 c??/'i (T (T m ). The Eq|311can be written as, 



S= [ drdr ^(r)(G° \ jmr >% e 
Jo j 

I J (7 (7 



{r) + S r 



(39) 



where, 



w*o> ( ^ - a + c' ) v w = jf dr ' x 



\ E(^^'^« _/x ) 



dr'(G° \ w (r-r')^ V (r') 



) 



(40) 
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and the clean non-interacting Green function matrix, 
G°(iid n ), is denned by 



The matrix element of Eq021 is given by the following 
Dysons like equation for the super-cell sites, 

G sc (I,J;ioj n )=G(I,J;iuj n ) 

L; iw n )E sc (i, L'; iuj n )G sc (L', J; iw n ). 



(41) 



L,L' 



By hint of Eq^| the clean Green function, G°, can be 
express in terms of the average Green function, G , and 
self energy, S, as, 

G 0_1 = G -1 + S. (42) 
By inserting Eq221 m to Eq[3H|we have, 

5 = | drdr J2 (G- 1 ) . w ^ v (r) + S mi 



(49) 



where, 

Sr — i — s 



(43) 



E J dr^oU™ (r)V^v W 

ijuu 
ijuu 

E J dr dr^Ur)^' ' (r - r')^ CT (r'). 



The second part of the right hand side of EqQ3 is the 
super-cell excluded effective medium partition function 
which is easily integrable, due to bilinearity of the Grass- 
mann variables. But the first part is the partition func- 
tion of the super-cell where is embedded in an effective 
medium environment, is not integrable directly due to 
four point Grassmann variable in its integrand. Note 
that for the case of only disordered systems, Eg 1471 leads 
to, 

G™ p ~\t - t) = g(iw n ) - 1 (r t') 5(t t') b sc (50) 

where e sc is the super-cell impurity matrix. The imagi- 
nary time Fourier transform of Eg 1501 imply that, 



(51) 



where the matrix element of Eo l51l can be written as, 
G% p (I,J,tu n ) = 9(I,J,tu n ) 



% j a a 



(44) 



+ E 9(J-, L, iu n )eLG l ™ p (L, J, uo n ) 

L 

(52) 



Now by applying the EMSCA 9 , where is taking average 
over all super-cells except one super-cell, which is denoted 
by {/}, Eg 1331 converts to, 

Zemsca =Z sc x JU miha (# i(7 #^) 

x e -J E„' E„«,, J} W V 

(45) 

where the super-cell partition function, Z sc , is given by, 

Z sc = (J nf = c f (d^ Ia d^, Ia ) e- s °-> ) r _ sc (46) 

and the super-cell action in the effective medium, S*^_n 
is, 

S S rU = E / 'drdr fair) (g-^jj^iT-T'WjAT) 

IJuu 

E / drM^U^ja' (r)ip J<T > (r) 

IJaa' 
IJcrcr 

in which the super-cell cavity Green function matrix, £?, 
is defined by, 



which are derived previously by us A Therefore, for a 
disordered system, Eos l32l49l andlCTare construct a close 
set of equations that should be solved self consistently. 

In the general case, where is included both interac- 
tion and randomness, to calculate the super-cell partition 
function ,Z SC , it is possible to use the Hirsch-Hubbard- 
Stratonovich transformation (HHST)A£ to decouple the 
interaction term and map it to an auxiliary Ising filed. 
Although our method is general with respect to the inter- 
action potential and randomness, but for simplicity and 
also to see the DCA derivation we concentrate our discus- 



sion on a repulsive on-site potential, U° 



and Stf? = 0. The HHST procedure is as following, di- 
viding the imaginary time interval, [0, /?], into M subin- 
tervals, At = hence the imaginary time at ith slice is 
given by 77 — Ijj. Therefore the discertizing of imaginary 
times leads to dr = J^i At—, thus, 

z sc = (f n CT nf = c ! I # Jct # Jct E ] e ~ s " )* 

V {«7=±l}/ 



where the super-cell action £*f^ is, 



(53) 



& 1 = G sc - S sc- 



(48) 



S?_i = (At) 2 E $iM >< 

IJll'a 
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(54) 

We define the cluster impurity Green function G™ p as 
(G*™ p = (G~ 1 )iju>+Sij{^^-^+^)5 ll > + 

The matrix form of Eg 1551 can be written as, 

(G™P- 1 )=g + T(e v "-l), (56) 
where, T JJH ' = ^ T yi and matrix elements of V CT are, 

V /i(T = Acts/; -e/Ar+ (57) 

Easl56landl55lfor the case of just interacting systems are 
derived by the DCA method^. Details of QMC solving 
of EqlSSl can be find in elsewher o 8 ! 15 . 

The average of impurity Green function, G*™ p , over 
Ising fields and also impurities configurations is the 
super-cell effective medium Green function 

(G i ™ p (I,J;T l ,T l ,)) = G sc (I,J;T h T f ). (58) 

The Fourier transform of Eg 1581 is give by, 

G sc (K n ;iu n ) = — x 

Yl e lK - r "e-"( T '-V)G sc (I, J; T U r r ). (59) 
11' u 

Eqs|2SI EH an d construct a close set of equations 
that should be solved self consistency. 

The algorithm of numerical proses is as following, 

1- A guess for the initial cluster self energy, S(i<T n ; iw), 
usually zero. 

2- From Eq|3Sl calculate the cluster Green function, 
G{K n] uj). 

3- Calculate the cavity Green function Q(K. n ;iuj) 
from Fourier transform of Eq0§| S _1 (K„; zu; n ) = 
G-^KnJMJn) + £(K„ ;«<;„). 



4- Calculate the Fourier transform of the cavity Green 
function, 

g(I,J;n-T r ) = i^^g(K„ ;i a;„)e^ e - iK " r " 

n K„ 

(60) 

5- calculate the new cluster Green function G(I, J; ri — 
Tf ) from EqEi ESI 

6- Calculate the inverse Fourier transform of 
G(J s J;7i-7-,/) 

7- Calculate the new self energy S(K n ; iu) n ) from 

S(K n ; iuj n ) = g-\K n - iw n ) - G- 1 ^; iu n ). (61) 

8- Go to step 2 and repeat whole proses until conver- 
gence. 

In conclusion we are extended the EMSCA to the 
case of an interacting disordered system. Similar to the 
EMSCA we showed that the periodicity of self energy 
with respect to super-cell translation vector leads to the 
coarse graining of self energies and hence the average 
Green function in k-space. Then by applying the effective 
medium theory on the system partition function, we find 
two equations where are relates the super-cell impurity 
Green function, G tmp sc (I, J; iui n ), and super cell average 
Green function, G sc (7, J; iui n ), to the super-cell cavity 
Green function, C/(J, J;iu) n )- This completes the whole 
formalism of a new real space method for disordered in- 
teraction systems. In especial case of interacting system 
our formalism leads to real space derivation of DCA 7 
while for a disoreded system is converts to the EMSCA 9 . 
Now we are established that the DCA and also the NL- 
CPA are two especial case of the, real space dynamical 
effective medium super-cell approximation (DEMSCA). 
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